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1 1 ABSTRACT 

O" 

We developed the tool GEM-FIND that allows to constrain the morphology and brightness distribution of ob- 
i— H ■ jects. The software fits geometrical models to spectrally dispersed interferometric visibility measurements in the 
| ^ j TV-band using the Levenberg-Marquardt minimization method. Each geometrical model describes the bright- 
ness distribution of the object in the Fourier space using a set of wavelength-independent and/or wavelength- 
dependent parameters. In this contribution we numerically analyze the stability of our nonlinear fitting approach 
by applying it to sets of synthetic visibilities with statistically applied errors, answering the following questions: 
How stable is the parameter determination with respect to (i) the number of itv-points, (ii) the distribution of 
points in the uw-plane, (iii) the noise level of the observations? 

CL[ 1. INTRODUCTION: GEM-FIND 

In order to fit observations obtained with the mid-IR interferometric instrument VLTI/MIDI with geometrical 
£j . models the software GEM-FIND (GEometrical Model Fitting for INterferometric Data) was developed. It fits 
^ 1 wavelength-dispersed visibility measurements (from 8-13 /im) with synthetic visibilities from centro-symmetric 
and asymmetric geometrical models. Each model represents the Fourier transformed brightness distribution of 
the object that is derived by means of a set of wavelength-independent and/or wavelength-dependent parameters. 
The wavelength-independent parameters are varied using an equidistant grid. For each grid point a non-linear 
least squares fitting minimization (based on the Levenberg-Marquardt method: see Ref.[l|) is performed for each 
. wavelength point on the wavelength-dependent parameters. This determines a color-reduced \ 2 using 

<n: x 2 

xl = ^f £ lor/ l...n,j 1....V. (1) 

o 

at grid point i and wavelength index j. Note that x\ j is already reduced for the number of degrees of freedom for 
• • ■ each wavelength point. The number of degrees of freedom in the fitting is solely determined by the wavelength- 
, ^ ' dependent parameters. The best-fitting parameters are then determined by 

Xmin = min (xi) ■ (2) 

& . 

Uncertainties of wavelength-independent parameters are derived from the 68.3% confidence region on each indi- 
vidual parameter using 

A x 2 =xl-xLn = A, (3) 
where A depends on the number of wavelength- independent parameters (see Ref.0). 

Up to now 13 models can be fitted with GEM-FIND. The different models and their parameters are given in 
Table HJ The flow-chart in Fig.Q] illustrates the fitting algorithm for GEM-FIND. 
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Table 1. Parametric description of the geometrical models included in GEM-FIND. 
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Notes. UD. .. uniform disk; CSE. . . circumstellar environment; d(maj) ■■■ uniform disk diameter (major axis); 
FWHM( mo:J ) . . . full width at half maximum of Gaussian (major axis); ^...inclination angle of ellipse; r\. .. axis ratio 
minor/major axis of ellipse; Aa;, Ay. . . angular offsets of binary component from center of symmetry; /. . . flux ratio bi- 
nary/primary or central star/envelope; 8 cen ... diameter of central star; 6 pllul , 8 SCC ... diameter of primary or secondary 
component; 9( ou t), 0(i n ). ■ ■ outer or inner ring diameter. 



2. MONTE-CARLO SIMULATIONS 

In order to test the reliability and limits of the geometrical models we derived synthetic visibilities with sta- 
tistically applied errors for a subsample of the models. These synthetic visibilities were then used as input 
quantities for GEM-FIND, i.e. we simulated a statistically significant number of observations (sample size 50) 
and performed a least-squares analysis for each of them. In this work we will present Monte-Carlo simulations 
for the models CircUD+CircGauss (M5), CircUD+EUGauss (M6) and UD+Dirac (M7) in TableQ] Models M5 
and M6 are optically thin and are especially suited to describe the circumstellar environment of e.g. asymptotic 
giant branch stars with low mass loss rate (e.g. Ref.0). Model M7, on the other hand, is able to describe any 
system composed of a resolved primary star and an unresolved secondary component. In the following we will 
describe the idea and goals of the Monte-Carlo simulations and the determination of the synthetic visibilities. 
Results are discussed in Sect.[3J 

2.1 Idea and Goals 

The Monte-Carlo simulations allow a numerical analysis of the stability of our fitting approach, thereby answering 
the following questions: 

• How stable is the parameter determination with respect to the number of uw-points N uv l 
The focus is on a scarce itv-coverage as the limited number of telescopes (2 in the case of MIDI) often only 
provides a sparse sampling of the ww-plane. We determined synthetic visibilities with N uv = 8, 6, 4 (upper 
row in Fig. [5]). 

• How does the distribution of w-points affect the result? 

We selected 4 different baseline configurations (middle row in Fig. [2]): (i) ideal (baseline lengths B p and 
position angles PA are uniformly distributed in the ww-plane), (ii) sameB (PA is uniformly distributed 
as in ideal, B p remains unchanged), (iii) samePA (B p is uniformly distributed as in ideal, PA remains 
unchanged), (iv) obs (we used the observed configuration from Ref.y). 
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Figure 1. Flow-chart illustrating the different stages in the fitting approach. Column A: The observed visibilities V 
are passed on to GEM-FIND and one of the geometrical models M1-M13 (see Table[JJ is selected. Column B: The 
wavelength-independent parameters ai are varied over an equidistant grid with n grid points. A Levenberg-Marquardt 
(LM) least squares fitting minimization is performed to the observed visibilities for each wavelength point Xj, i.e. 



the 
A 



wavelenght-dependent parameters bij are optimized to minimize \ ■ The best solution results in a reduced Xi,j- 
color-reduced xf is derived for each grid point. Column C: The best-fitting wavelength-dependent and independent 
parameters b' and a' are determined by the grid point i having the minimal Xmin- 



• How does observational noise influence the result? 

Different noise levels at FWHM a were applied to the synthetic visibilities. As shown in Ref.0 the typ- 
ical noise level of MIDI observations lies between 7% (very good quality) and 15% (acceptable quality). 
Therefore, we statistically applied Gaussian distributed errors of a = 7, 10, 15% to the synthetic visibilities 
(lower row in Fig. [5]). 

• How are the errors of the resulting parameters distributed? A least squares fitting minimization 
provides not only the best fitting values of the parameters but also the covariance matrix of the estimates 
that allows to calculate the uncertainties on the fitted parameters. These uncertainties are reliable if the 
errors on the parameters are known to be normally distributed. In order to guarantee this, observational 
errors have to be Gaussian distributed and the model must be linear. In certain cases even for a nonlinear 
model uncertainties can be derived from the covariance matrix (Ref.0), but this has to be tested using 
a Monte-Carlo simulation. If uncertainties are not Gaussian distributed, appropriate confidence regions 
have to be constructed. We determine the error of the results by calculating the difference between the 
true parameter value and the value calculated with GEM-FIND (summed over wavelength for wavelength- 
dependent parameters). To judge the level of normality a \ 2 test is applied. This test provides the 
probability of getting the calculated x 2 for Gaussian distributed errors. Our null hypothesis is that errors 
are Gaussian distributed, i.e. a x 2 probability of 0.05 (this corresponds to a typical significance level of 
5%) or lower points to a statistically significant deviation from this hypothesis. The x 2 probability is given 
in the upper right corner in the insets of Fig.l4lfT2l 

Details on the synthetic visibility sets are given in Tabled 
2.2 Synthetic visibility determination 

With interferometry we observe the complex visibility of an object, i.e. the Fourier transform of its brightness 
distribution. The Earth's atmosphere influences the observations and makes it difficult to access the absolute 
phase of the object. Therefore, in the following only the absolute value of the complex visibility, the normalized 
visibility, will be considered. Synthetic visibility functions can be derived from analytical formulae in the uv- 

B B 

plane, with u = -f - sin (PA) and v = -f - cos (PA). The synthetic visibility of model M5 can be derived by 
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Figure 2. Upper row: iro-coverage with N uv = 8, 6, 4 (from left to right) used to determine synthetic visibilities. 
Contour levels range from 8-12.5 /im (light to dark, respectively) with a step size of 0.5 /im. Middle row: it«-coverage 
of the configurations ideal, sameB, samePA, obs (from left to right) used to determine synthetic visibilities. Lower row: 
Examples of synthetic visibilities versus wavelength with different noise levels added (7, 10, 15%; from left to right). 



with Ji being the first order Bessel function of the first kind. The value of r = yju 2 + v 2 depends on the 
configuration (ideal, sameB, samePA, obs) that is used. 

Applying a rotation and compression to one axis (which becomes the semi-minor axis) of a circular model results 
in an ellipse, i.e. 

u,/, = u cos ip — v sin ip and v^, = u sin ip + v cos ip , (5) 



Tip, 7] 



(6) 



where ip denotes the inclination angle of the ellipse major axis (in sky plane, or minor axis in Fourier plane) 
measured from North to East and rj is the ratio of minor over major axis dicim6t6r tj — @min /9maj ■ The visibility 
function for model M6 is then given by 
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The synthetic visibility of model M7 can be derived with 

2Jl J r T 9prim) + / exp^TTz^Ax + vAy)) 
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1 + / 



(8) 



Table 2. Details on the synthetic visibility sets. Given is the reference number of the synthetic visibility set, the number 
of Mii-points, the type of configuration and the assumed noise level. 
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A wavelength- independent, random relative error oy is statistically applied to the synthetic visibilities by means 
of a random number generator and results in 

Vr(u,v)=Vi(u,v)+a r , (9) 

with i defining the model that is used. The probability distribution of the error is assumed to be normal with a 
FWHM of either 7%, 10% or 15%. 

The selection of the input parameters for Vi will be discussed in the following. The parameters below are set in 
a way to describe a typical environment of an AGB star. The geometrical models implemented in GEM-FIND 
can be easily applied to other objects as well (e.g. YSO, AGNs). 

2.2.1 Choice of wavelength-independent parameters 

CircUD+CircGauss: The diameter of the central star is fixed to 8 cen = 10 mas. The diameter of the central star 
is in practice often assumed to be known and therefore enters as a fixed parameter. 

CircUD+EllGauss: The diameter of the central star is fixed to 9 cen = 10 mas, the axis ratio rj = 0.5 and the 
position angle of the ellipse ip = 45°. 

UD+Dirac: Observed separations of binary systems involving an AGB star lie between 7 mas and 5" (see Ref.Q). 
In this work the separation s of the companion (i.e. s = y Ax 2 + Ay 2 ) is fixed to s = 40 mas. 

2.2.2 Choice of wavelength-dependent parameters 

CircUD+CircGauss: The wavelength-dependent flux ratio / as well as the FWHM of the Gaussian envelope 
were determined with the model of RR Aql (Fig. 10 of Ref. 7), a well studied oxygen- rich Mira variable. Ref.@ 
studied the pulsation mechanism of the star by modeling multi-epoch interferometric observations with a radiative 
transfer code. We set the flux ratio by dividing the flux of the central star by the flux of the envelope (upper 
panel in figure). The FWHM was set by fitting a CircUD+CircGauss model to the visibilities (lower panel in 
figure). The determined parameters are plotted in the upper row of Fig. [3] 
CircUD+EllGauss: Parameters were set in the same way as described above. 

UD+Dirac: The diameter of the uniform disk was set as described above, but instead of a CircUD+CircGauss 
model a UD model was fitted to the visibilities (lower panel in Fig. 10 of Ref.0). The flux ratio between the 
binary companion and the primary was set with the help of two blackbodies. One blackbody with a temperature 
of 400 K (mean temperature of dust shell, Ref.@) was used to set the flux of the primary. A secondary component 
(e.g. planet) that is present in the surrounding of an AGB star will accrete matter from the mass that is lost by 
the star. Therefore, the other blackbody with a temperature of 300 K (depending on the separation, Ref.|9|) was 
used to set the flux of the companion. Determined parameters are shown in lower row of Fig. [3] 
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Figure 3. Upper row: Parameters used for the synthetic visibility determination of CircUD+EUGauss and Cir- 
cUD+CircGauss. Left: Wavelength-independent axis ratio rj for CircUD+EUGauss (full line) and CircUD+CircGauss 
(dashed line). Middle: Wavelength-dependent flux ratio /. Right: Wavelength-dependent FWHM of the envelope (dark 
violet line) as well as fixed diameter of central star (light blue line). Lower row: Parameters used for the synthetic 
visibility determination of UD+Dirac. Left: Wavelength-independent position of the unresolved companion. Middle: 
Flux ratio / of the companion over the primary. Right: Wavelength-dependent diameter of the primary. 

3. RESULTS 

The derived sets of synthetic visibilities (S1-S36, Table© were used as input for GEM-FIND for models M5, M6 
and M7. In the following we describe the result of the fit for each model separately, answering the questions 
raised in Sect. 12. II 

3.1 CircUD+CircGauss 

The CircUD+CircGauss model is spherically symmetric and has 3 input parameters (flux ratio /, FWHM of the 
Gaussian, central star diameter 9 C en), where the first two are wavelength-dependent. The diameter of the central 
star is often assumed to be known and is therefore considered as a fixed wavelength-independent parameter, i.e. 
the effective number of free parameters (i is 2. 

How stable is the parameter determination with respect to the number of uu-points N uv l Fig.|Hplots the results 
obtained with GEM-FIND for the synthetic visibility sets SI, S13, S25 (from left to right): FWHM (upper pan- 
els) and flux ratio / (lower panels) versus wavelength. The insets show the error distribution in the sample 
between the true value and the values derived with GEM-FIND. The results for the flux ratio as well as the 
FWHM are comparable to the true value for N uv = 4, 6, 8. One would expect that the error distribution on the 
parameters is narrower the more uv-points are used. But, as can be seen from the figure, the RMS for 6 points 
is larger than the one for 4 points. This behavior is caused by the difference in the baseline lengths between the 
4 points configuration (B p =120, 120, 40, 40 m) and the 6 points (B p =120, 120, 70, 70, 20, 20 m) configuration. 
The smallest baseline in the 6 points configuration is resolving the star only slightly. The visibility of these 
points is very close to 1 and therefore these points do not add much information, i.e. we are left with 4 points 
also for the 6 points model. As the separation of the remaining baseline lengths is smaller than for the 4 points 
model, the determination of the right FWHM is more difficult. 

To summarize, for this spherically symmetric model the parameter determination is stable for N uv — 8,6,4. 

How does the distribution of uu-points affect the result? Fig. [5] plots the parameters obtained with GEM-FIND 
for the synthetic visibility sets S25, S28, S31 and S34 (from left to right). All configurations except sameB 
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Figure 4. Parameters obtained with GEM-FIND for the synthetic visibility sets SI, S13, S25 (from left to right) of model 
CircUD+CircGauss: FWHM (upper panels) and flux ratio / (lower panels) versus wavelength. The blue line corresponds 
to the true value. Wavelength-averaged means and standard deviations are denoted. The insets show the error distribution 
in the sample. The probability of normality is given in the upper right corner. 



(second column) find the right solution. sameB uses the same baseline length for all points, i.e. in a spherical 
environment, all the points carry the same information (up to the synthetic errors). Therefore, the wu-coverage 
is too small to find the right solution. The similarity between the results of configurations ideal (first column) 
and samePA (third column) is expected for a spherically symmetric model, because the visibility profile looks 
the same no matter which position angle is used. The configuration obs (fourth column) has a larger RMS for 
the flux ratio /, but is almost the same for the FWHM. This effect can be explained by the larger range in 
baseline length for configurations ideal and samePA, which is of greater importance for the determination of the 
flux ratio /. 

In summary, the ideal and samePA are most suitable to find the right parameters for the spherically symmetric 
CircUD+CircGauss model. Also the obs configuration is able to determine the correct parameters, but with a 
slightly broader error distribution. The sameB configuration is not usable for spherically symmetric models with 
more than one component. 

How does observational noise influence the result? Fig. [5] plots the parameters obtained with GEM- FIND for the 
synthetic visibility sets S25, S26, S27 (from left to right). As expected, the larger the noise level on the visibilities 
the larger the RMS. The right solution is also found for a noise level as high as 15%, even if the FWHM and 
flux ratio / tend to be overestimated/underestimated, respectively. An even higher uncertainty may already be 
critical. 

How are the errors of the resulting parameters distributed? The insets in Fig. 21 Fig. [5] and Fig. H2 show the error 
distribution in the sample (normalized and summed over all wavelengths) between the true value and the value 
calculated with GEM-FIND. The probability of normality is given in the upper right corner. The wavelength- 
averaged means and standard deviations given in the figures are very close to the true value FWHM tr uc=28.0 
and /true=0.60. The error distributions for the flux ratio are approximately Gaussian (probability >0.05) for 
all of the configurations except sameB. On the other hand, error distributions for FWHM are not Gaussian 
distributed (probability <0.05). The FWHM tends to be slightly overestimated. 



3.2 CircUD+EllGauss 

The asymmetric CircUD+EllGauss model has 5 input parameters (flux ratio /, FWHM of the Gaussian, central 
star diameter 8 cen , inclination angle of the ellipse ip, axis ratio 77), where the first two are wavelength-dependent. 
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Figure 5. Same as Fig. [4] but for synthetic visibility sets S25, S28, S31 and S34 (from left to right) of model Cir- 
cUD+CircGauss. 
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Figure 6. Same as Fig. [4] but for synthetic visibility sets S25, S26, S27 (from left to right) of model CircUD+CircGauss. 



100 

80 

60 

40 

20 


3.5 

3.0 
2.5 
2.0 
1.5 
1.0 
0.5 
0.0 

1.0 
0.8 
0.6 
0.4 

0.2 

0.0 



4points ideal 7% 



6points ideal 7% 



8points ideal 7% 




Mean 0.58 
RMS 0.05 




Mean 28.3 
RMS 0.6 




-4 -2 2 4 
Error [a] 




10 



1 1 



12 



9 10 11 

Wavelength [/Am] 



12 8 



10 



12 



30 r 

Meon_V 44.3 25 
RMS-V 7-2 20 
" Meon_Tj 0.49 z 15 
RMS_T) 0.02 '0 



-4 -2 2 4 
Error [a] 



Mean_V 44.6 
RMS_V 4 - 7 
Meon_Tj 0.49 
RMS_7j 0.04 

+ + 



30 r 

Meon_V 44.8 25 - 

RMS-V 1-6 20 - 

" |- Mean_7j 0.50 z '5 
RMS_jj 0.00 1 ° " 



-4 -2 2 4 
Error [a] 



-4 -2 2 4 
Error [u] 



o 



50 



100 



150 



50 



100 



150 



50 



100 



150 



t 



Figure 7. Parameters obtained with GEM-FIND for the synthetic visibility sets SI, S13, S25 (from left to right) of model 
CircUD+EUGauss. Upper+middle: FWHM (upper panels) and flux ratio / (middle panels) versus wavelength. Lower: 
Axis ratio r\ versus inclination angle ijj. The blue line/cross corresponds to the true value. Means and standard deviations 
are denoted. The insets show the error distribution in the sample. The probability of normality is given in the upper 
right corner. 



The diameter of the central star enters as a fixed parameter. 

How stable is the parameter determination with respect to the number of uu-points N uv ? Fig.[7]plots the param- 
eters obtained with GEM-FIND for the synthetic visibility sets SI, S13, S25. For the FWHM (upper panels) and 
flux ratio / (middle panels) the values found with GEM-FIND are very close to the true values for N uv = 4, 6, 8. 
As anticipated, the RMS is larger the less points are used. The observed behavior for 6 points described in 13.11 
appears also in the case of an ellipse (evident from the mean and standard deviation values given in the figure) . 
The correct solution for the axis ratio rj (y-axis in lower panels) is found for N uv = 4, 6, 8. On the other hand, 
the correct value for the inclination angle ip of the disk is only found in case of N uv = 8 and to a lesser extent 
also in case of N uv = 6. 

How does the distribution of uw-points affect the result? Fig. [5] plots the parameters obtained with GEM-FIND 
for the synthetic visibility sets S25, S28, S31 and S34. The ideal (first column) and obs (last column) configu- 
rations are both able to reliably determine the right value for all parameters. The configuration samePA (third 
column) is not able to retrieve FWHM, r\ and i/j, because using the same position angle for all points does not 
result in enough independent data points to determine the parameters of an ellipse. On the other hand, the 
flux ratio / does not depend on the ellipticity (and consequently not on the position angle) and can therefore 
be determined by this configuration. The configuration sameB (second column), on the other hand, is able to 
determine FWHM, r/ and ip. The flux ratio / is more difficult to determine, because differences in flux ratio cause 
a more pronounced difference in the visibility profile at larger baselines. The larger baseline lengths (~12Q m) 
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Figure 8. Same as Fig. [7] but for synthetic visibility sets S25, S28, S31 and S34 (from left to right) of model Cir- 
cUD+EHGauss. 



are not included in the configuration sameB (largest baseline ~60m), but covered by ideal and samePA. The 
largest baselines for obs are ~90m which accounts for the slightly larger RMS. 

How does observational noise influence the result? Fig. [3] plots the parameters obtained with GEM-FIND for the 
synthetic visibility sets S25, S26, S27. The correct parameters are detected for 7%, 10% as well as 15%. As 
expected, the larger the noise level on the synthetic visibilities the larger the RMS. 

How are the errors of the resulting parameters distributed? The error distribution is most of the time non-Gaussian 
for all parameters (probability <0.05). This suggests that we cannot use the errors provided by the covariance 
matrix. Due to the coarse step size that is used for the axis ratio r\ (A77 = 0.1) it is not possible to determine the 
level of normality of the error distribution. Appropriate confidence regions should be constructed to determine 
the errors on all parameters. 

3.3 UD+Dirac 

The UD+Dirac model has 4 input parameters (position of the binary component Ax and Ay, flux ratio /, 
diameter of the primary # pr i m ), where the latter two are wavelength-dependent. 

How stable is the parameter determination with respect to the number of uu-points N uv ? Fig. [10] plots the pa- 
rameters obtained with GEM-FIND for the synthetic visibility sets S4, S16, S28. The correct values for the 
diameter and position of the companion can be found when 6 or more points are used. With 4 points the param- 
eter determination is no longer reliable (this is evident from the large RMS). Note that the flux ratio can also 
be determined with 6 and 8 points, but with an uncertainty as large as ±25%. The flux ratio in our samples is 
quite small leading to a low contrast in visibilities. A larger flux ratio likely would result in a lower uncertainty. 
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Figure 9. Same as Fig.[7]but for synthetic visibility sets S25, S26, S27 (from left to right) of model CircUD+EUGauss. 



How does the distribution of uu-points affect the result? Fig.[TT]plots the parameters obtained with GEM-FIND 
for the synthetic visibility sets S25, S28, S31 and S34. The ideal (first column) and samePA (third column) 
configurations are not able to determine the expected values for the diameter, flux ratio and position of the com- 
panion. The reason for this behavior is the different baseline lengths, i.e. for the binary model it is advantegous 
to use baseline lengths as similar as possible. This statement is confirmed by the sameB configuration and also 
by the obs configuration that show only a small variation in baseline lengths. These configurations are able to 
find the right solution for all parameters. 

How does observational noise influence the result? Fig.[T2] plots the parameters obtained with GEM-FIND for 
the synthetic visibility sets S28, S29, S30. Noise levels of 7% and 10% are able to determine the right param- 
eters. Also a noise level of 15% is still able to determine the diameter and position of the companion with a 
reasonable error distribution. On the other hand, the determination of the flux ratio is no longer possible with 
such a high noise level. As mentioned before, we are confident that a larger value of the flux ratio would allow 
to find the right solution with a lower uncertainty. 

How are the errors of the resulting parameters distributed? The error distribution of the parameters is non- 
Gaussian for most synthetic visibility sets. Also in this case it is recommended to use confidence regions for the 
error estimation. 
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Figure 10. Parameters obtained with GEM-FIND for the synthetic visibility sets S4, S16, S28 (from left to right) of 
model UD+Dirac. Upper+middle: UD diameter (upper panels) and flux ratio / (middle panels) versus wavelength. 
Lower: x-position versus y-position of spot. The blue line/cross corresponds to the true value. The insets show the error 
distribution in the sample. The probability of normality is given in the upper right corner. 



4. CONCLUSIONS 

In this work we presented the geometrical model fitting tool GEM-FIND for MIDI interfcromctric data. In order 
to guarantee reliability and to test the limits of the fitting method a Monte-Carlo approach was used to generate 
synthetic visibilities with statistical noise of three models (CircUD+CircGauss, CircUD+EllGauss, UD+Dirac). 
We used synthetic data sets with different itu-coverage, number of uw-points and noise levels in order to be able 
to answer the following questions: How many uu-points are needed for a correct parameter determination? How 
do these points have to be distributed in the iw-plane? How large can the noise level of the observations be? 
The following main results were found: As expected, it is recommended to (i) avoid the use of similar baseline 
lenghts (sameB) to determine the parameters of a spherical object morphology; (ii) avoid the use of similar 
position angles (samePA) to determine the parameters of an elliptical object morphology, (iii) We find that the 
use of close baseline lengths (sameB and obs) facilitates the determination of the parameters of a binary system. 
Here great care must be taken when the morphology of an object with unknown geometry is to be determined, 
i.e. one has to account for the effect of distribution of uw-points on the fitting results, (iv) Our study shows that 
observations with a noise level as high is 15% can still be used to determine the parameters in most of the cases, 
(v) As the error distribution of the output deviates from normality, errors of all parameters should be derived 
using confidence regions. 
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Figure 11. Same as Fig.QIJ]but for synthetic visibility sets S25, S28, S31 and S34 (from left to right) of model UD+Dirac. 
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